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Abstract 

The question about the existence and characterization of bounded solutions 
to linear functional-differential equations with both advanced and delayed ar- 
guments was posed in early 1970s by T. Kato in connection with the analysis of 
the pantograph equation, y'(x) = ay(qx) + by(x). In the present paper, we an- 
swer this question for the balanced generalized pantograph equation of the form 
-a 2 y"(x) + aiy'(x) + y(x) = J* °° y(ax) fJ.(da), where ai > 0, a 2 > 0, a\ + a\ > 0, 
and n is a probability measure. Namely, setting K := / °° In a //(da), we prove 
that if K < then the equation does not have nontrivial (i.e., nonconstant) 
bounded solutions, while if K > then such a solution exists. The result 
in the critical case, K = 0, settles a long-standing problem. The proof ex- 
ploits the link with the theory of Markov processes, in that any solution of 
the balanced pantograph equation is an ^-harmonic function relative to the 
generator £ of a certain diffusion process with "multiplication" jumps. The pa- 
per also includes three "elementary" proofs for the simple prototype equation 
y'(x) + y(x) = \ y{qx) + \ y{x/q), based on perturbation, analytical, and prob- 
abilistic techniques, respectively, which may appear useful in other situations 
as efficient exploratory tools. 

Key words: Pantograph equation, functional-differential equations, integro- 
differential equations, balance condition, bounded solutions, WKB expansion, 
(/-difference equations, ruin problem, Markov processes, jump diffusions, in- 
harmonic functions, martingales. 
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1. Introduction 



The classical pantograph equation is the linear first-order functional-differential equa- 
tion (with rescaled argument) of the form 

y'(x) = ay(qx) + by(x), (1) 

where a, b are constant coefficients (real or complex) and q > is a rescaling parame- 
ter. Historically^] the term "pantograph" dates back to the seminal paper of 1971 by 
Ockendon and Tayler [23] , where such equation^] emerged in a mathematical model 
for the dynamics of an overhead current collection system on an electric locomotive 
(with the physically relevant value q < 1). At about the same time, a systematic 
analysis of solutions to the pantograph equation was started by Fox et al. [H] , where 
various analytical, perturbation, and numerical techniques were discussed at length 
(for both q < 1 and q > 1). 

Interestingly, an equation of the form (CQ) (with q > 1) was derived more than 25 
years earlier by Ambartsumian [2] to describe the absorption of light by the interstellar 
matter. Some particular cases of Eq. (pQ) are also found in early work by Mahler 
[19] on a certain partition problem in number theory (where Eq. (Ill) , with a — 1, 
b = 0, q < 1, appears as a limit of a similar functional-difference equation) and by 
Gaver [T2] on a special ruin problem (with a = 1, b = — 1, g > 1). Subsequently, 
the pantograph equation has appeared in numerous applications ranging from the 
problem of coherent states in quantum theory [26J to cell-growth modeling in biology 
[28] (see further references in Refs. [31 [71 [HI [20]). These and other examples suggest 
that, typically, the pantograph equation and similar functional-differential equations 
with rescaling are relevant as long as the systems in question possess some kind of 
self-similarity. 

Since its introduction into the mathematical literature in the early 1970s, the 
theory of the pantograph equation (and some of its natural generalizations) has been 
the subject of persistent attention and research effort, yielding over years a number 
of significant developments. In particular, the classification of Eq. ([T]) with regard 
to various domains of the parametersjl including existence and uniqueness theorems, 
and an extensive asymptotic analysis of the corresponding solutions have been given 
by Kato and McLeod [T7j and Kato [16] . The investigation of such equations in the 
complex domain was initiated by Morris et al. [22] and Oberg [23] and continued 
by Derfel and Iserles [H] and Marshall et al. [20]. A systematic treatment of the 
generalized first-order pantograph equation (with matrix coefficients and also allowing 
for a term with rescaled derivative) is contained in the influential paper by Iserles [13] , 
where in particular a fine geometric structure of almost-periodic solutions has been 
described. Asymptotics for equations with variable coefficients have been studied by 
Derfel and Vogl [5] ■ 

lr rhe name "pantograph equation" was not in wide use until it was coined by Iserles |14) for a 
more general class of functional-differential equations. 
2 To be more precise, a certain vector analog of Eq. |TJ). 

3 Depending on whether q < 1 or q > 1 and also on the cases 5R6 < 0, 5R6 > 0, and = 0. 
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Higher-order generalizations of the pantograph equation (CQ) lead to the class of 
linear functional-differential equations with rescaling, 



EEW*W + &)=0 (2) 

j=l k=0 

(see Ref. [5] and further references therein)!^ Kato [16j posed a problem of asymptotic 
analysis of Eq. (J2J), including the question of existence and characterization of bounded 
solutions. Some partial answers to the latter question have been given by Derfel |U [5] 
and Derfel and Molchanov [7|. 

In particular, Derfel |1] considered the "balanced" generalized first-order panto- 
graph equation of the form 

e 

y'(x) +y(x) = ^pjyiajx), (3) 

3=1 



subject to the condition 



^2 Pi 



Pj >0 {j = !,...,£), (4) 



so that the weights P j of the rescaled |/-terms on the right-hand side of Eq. ([3]) match 
the unit coefficient of the y(x) on the left. Note that, owing to the balance condition 
fllD, Eq. ([3]) always has a trivial solution y = const. The question of existence of 
nontrivial (i.e., nonconstant) bounded solutions is most interesting (and most difficult) 
in the case where the right-hand side of Eq. involves both "advanced" (ay > 1) 
and "delayed" (0 < a 3 - < 1) arguments. It turns out that the answer depends crucially 
on the quantity 

t 

K ^J^Pjlnaj. (5) 

3=1 

Namely, Derfel [I] has proved that if K < then Eq. ([3]) has no nontrivial bounded 
solutions, whereas if K > then such a solution always exists. In the "critical" case 
K = 0, this question has remained open as yet. 

In the present paper, we consider a more general integro-differential equation^! of 
the pantograph type, namely, 

POO 

- a 2 y"(x) + aiy'(x) + y(x) = / y(ax) fi(da), (6) 

J o 



4 Note that theory of such equations is closely related to the theory of g-difference equations 
developed by Birkhoff [3j and Adams [1] (see also Section [3] below) . 

5 In fact, the results of the paper [4] mentioned above include first-order equations of the form 
(|6j). i.e., with 02 = 0. Let us also remark that more general first-order integro-differential equations 
(but with delayed arguments only, i.e., a S (0, 1)) were considered by Iserles and Liu [15] . 
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where a± > 0, a 2 > 0, af + > (so that ai, 02 do not vanish simultaneously), and 
li is a probability measure on (0, 00), 

/■oo 

/i(0, 00) = / /u(da) = 1. (7) 



The parameter a in Eq. ([6]) can be viewed as a random variable, with values in (0, 00) 
and the probability distribution given by the measure /z, i.e., P{a G A} = fi(A), 
A C (0,oo). Note that Eq. is balanced in the same sense as Eq. ([3]), since the 
mean contribution of the distributed rescaled term y(ax) is matched by that of y(x). 
Moreover, Eq. §6§ reduces to Eq. ([3]) when ai = 1, a 2 = 0, and the measure /1 is 
discrete, with atoms pj = f-i(atj), j = !,...,£ (i.e., a is a discrete random variable, 
with the distribution P{a = ay} = Pj, j = !,...,£). 

As already mentioned, due to the balance condition ([7]) any constant satisfies 
Eq. (Q, and by linearity of the equation one can assume, without loss of generality, 
that y(0) = 0. Moreover, if x > (x < 0) then the right-hand side of Eq. (JH]) 
is determined solely by the values of the function y(u) with u > (respectively, 
u < 0). Therefore, the two-sided equation §6§ is decoupled at x = into two one- 
sided boundary value problems, 



-a 2 y"{x) + axy^x) + y(x) = / y(ax) fi(da), x ^ 0, ^ 







2/(0) = 0. 

For Eq. §6§ , the analog of Eq. §5§ is given by 

POO 

K := In a //(da) = E [In a). (9) 



Our main result is the following theorem, which resolves the problem of nontrivial 
bounded solutions in the critical case, K = (and also recovers and extends the 
result of Ref. [1] for the case K < 0). 

Theorem 1. Assume that 7^ E|lna| < 00, so that K in Eq. (TJ|) is well defined 
and the measure fi is not concentrated at the point a = 1, i.e., the random variable a 
does not degenerate to the constant 1. Under these hypotheses, the condition K < 
implies that any bounded solution of equation (0|) is trivial, i.e., y(x) = const, 

The apparent probabilistic structure of Eq. ([6]) is crucial for our proof of this 
result. The main idea is to construct a certain diffusion process X t , with negative 
drift and "multiplication" jumps (i.e., of the form x 1— > ax), such that Eq. ([6]) can 
be rewritten as Cy = 0, where £ is the infinitesimal generator of the Markov process 
X t . That is to say, the class of bounded solutions of Eq. coincides with the set of 
bounded £-harmonic functions. This link brings in the powerful tool kit of Markov 
processes; particularly instrumental is the well-known fact (see, e.g., Ref. [9]) that for 
any £-harmonic function f(x), the random process f(X t ) is a martingale, and hence, 
for any t > 0, 

f(x) = E[f(X t )\X = x], xeR. (10) 



4 



On the other hand, due to the multiplication structure of independent consecutive 
jumps of the process X t , its position after n jumps is expressed in terms of a back- 
ground random walk Sk = £1 + • • • + £jfc (0 — ^ — n )> where £»'s are independent 
random variables with the same distribution as In a. The hypothesis K < of The- 
orem [1] implies that, almost surely (a.s.), the random walk S n travels arbitrarily far 
to the left. Using an optional stopping theorem (whereby the boundedness of f(x) 
is important), we can apply the martingale identity (flO]) at the suitably chosen stop- 
ping (first-passage) times, which eventually leads to the conclusion that f'(x) = 
and hence f(x) = const. 

This approach also allows us to give an example of a nontrivial bounded solution 
to equation ([6]) in the case K > (thus extending the result by Derfel jl] to the 
second-order pantograph equation). 

Theorem 2. Suppose that K > 0, and set 

/oo(x) := Pjliminf X t = +oo\X = x\, x G E, (11) 

where X t is the random process constructed in the proof of Theorem Ql Then the 
function foa(x) is L-harmonic and such that foo(x) — > as x — > — oo and foo{x) — > 1 
as x — >• +oo. 

In the particular case where 02 = 0, Eq. ([6]) becomes 

POO 

aiy'(x) + y(x) = / y(ax) /i(da), 16R. (12) 
Jo 

In this case, the diffusion component of the random process X t is switched off, and it 
follows, due to the negative drift and multiplication jumps (see details in Section @J, 
that if X = x < 0, then X t < for all t > 0. That is to say, the negative semi-axis 
(—00, 0] is an absorbing set for the process X t , and hence the function foo(x), defined 
by Eq. (JTTJ as the probability to escape to +00 starting from x, vanishes for all x < 0. 
This leads to the following interesting specification of the example in Theorem [2j 

Corollary. If a 2 = in Eq. then foo{x) = for all x < 0. Moreover, Eq. 
(0|) implies that all derivatives of the function foo(x) vanish at zero, f^(0) = 
(fc = l,2,...). 

Before elaborating the ideas outlined above, we would like to make a short digres- 
sion in order to consider the simple prototype example of Eq. ([3]) , namely, 

1 1 

y'(x) +y(x) = -y(qx) + -y(q- l x) (q ^ 1), (13) 

and to give several different "sketch" proofs of Theorem [1] in this case. Note that, 
according to Eq. ([5]), we have 

K = - In q + - In - = 0, 
2 H 2 q 
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so Eq. ( Ti~3l) falls in the (most interesting) critical case. In fact, this example was 
the starting point of our work and a kind of mathematical test-tube to try various 
approaches and ideas. Although not strictly necessary for the exposition, after some 
deliberation we have cautiously decided to include our early proofs (based on per- 
turbation, analytical, and probabilistic argumentsj§ respectively), partly because this 
will hopefully equip the reader with some insight into validity of the result, and also 
because these methods may appear useful as exploratory tools in other situations. 

The rest of the paper is laid out as follows. In Sections [2J EH and@], we discuss the 
three approaches to equation ffl3|) as just mentioned. In Section [5] we start a more 
systematic treatment by describing the construction of a suitable diffusion process 
with multiplication jumps. In Section [61 we discuss the corresponding £-harmonic 
functions and obtain an a priori bound for the derivative of a solution. Finally, in 
Section [7] we prove our main Theorems [T] and [2J 

2. Perturbative proof 

Following the ideas used by Ockendon and Tayler [21] and Fox et al. [11] in the case of 
the original pantograph equation ([1]), we start by observing that if q — 1 then Eq. ffl3|) 
is reduced to the equation y' = 0, which has constant solutions only. Therefore, when 
the parameter q is close to 1, it is reasonable to seek solutions of Eq. (fT5|) in a form that 
involves "superposition" of (exponentially) small oscillations (fast variation) on top 
of an almost constant (polynomial) function (slow variation). This leads to a WKB- 
type asymptotic expansion of the solution in terms of perturbation parameter e ~ 
(see Ref. [2 5 J ) , which in the first-order approximation yields two first-order differential 
equations: a nonlinear equation (called the eikonal equation) for the fast variation and 
a linear equation (called the transport equation) for the slow variation. For simplicity 
of presentation, we will restrict ourselves to the first-order approximation, but in 
principle one can go on to the analysis of higher-order terms, which are described by 
linear equations and therefore can be determined without much trouble. 

To implement this approach, set q = 1 ± e (e > 0), x = e~ x u, and y(x) = 
y(e^ 1 u) =: f(u). Then 



As explained in the Introduction, without loss of generality we may assume that 



6 It is amusing that these three methods represent nicely the traditional organization of British 
mathematics into Applied Mathematics, Pure Mathematics, and Statistics, which is reflected in the 
names of mathematical departments in most universities in the UK. 




and Eq 



f|T3|) takes the form 



ef(u) + f(u) = l - /((l ± e)u) + l - / ((l ± e)-V). 



(14) 



/(0) = 0. 
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Now, suppose that, for small e, the function /(it) admits a WKB-type expansion, 

f(u) ~ (A (u) + eA 1 (u) + ■■■) exp( £ -V(«)), (15) 

which in principle should be valid uniformly for all u, including the limiting values 
u — > and u — ■> oo. Differentiation of Eq. ( fl5i) yields 

e/'(u) ~ (Aq(u)V'(u) +e(A' (u) + A^V'iu)) + • • • ) exp(^V( M )). (16) 

From Eq. ( TT5i) we also obtain 

/((l±e)«) ~ (A (u)+e(A 1 (u)±uA' (u)) + ---) 

fV(u) , T/// . aV» \ (17) 
x exp f — ± uV'(ii) H 1 1 

and 

f((l±e)- 1 u) ~ (A ( M )+ £ (A 1 ( M ) TM A[ ) ( M )) + ---) 

x exp — ^ =f uV(u) + £ uV'(u) H 7 x 1 



Substituting the expansions (flB]) . (iT7j) and (|T8|) into Eq. fTl4|) . canceling out the com- 
mon factor exp(V(u)/£), and collecting the terms that remain after setting e = 0, we 
get 

A (u)V'(u) +A Q (u) = ~A (u)exp(±uV'(u)) + ^A (u)exp(TuV'(u)). 

Assuming that Aq(u) ^ 0, this gives the equation 

1 + V'(u) = cosh(uV'(u)), (19) 

or equivalently 

uV'iu) . . 

u = -. Ht . 20 

cosh(uV'(u)) - 1 V ; 

Similarly, equating the terms of order of e and noting that Ai(u) cancels out owing 
to Eq. ( fl9l) . we obtain 

A' (u) = A' {u) u siiih(u7'(u)) + - A (u) u 2 V"{u) cosh{uV {u)) (21) 

+ - A Q (u) uV'(u) exp(=F«\/'(w)) 

We can now check that the formal expansion (|15|) is compatible with the zero 
initial condition, /(0) = 0. Equation (j20|) implies that if u — > then uV'(u) — > oo, 
and moreover 

uV'(u) ~ -lnw + lnln-H , (22) 
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whence 



ln 2 w 1 
V{u) ~ — + lnu-lnln~ + --- . (23) 



Furthermore, differentiation of Eq. f|22l) gives 

u 2 V"(u) ~lnu-lnln- + --- . (24) 

Inserting formulas (122]) and (!24l) into Eq. (I2T1) . we obtain for Ao(w) the asymptotic 
differential equation 

~ — In u - In In - + • • 



Aq{u) 2u V u 
which solves to 

In it / 1 \ 

\nA (u) flnw-21nln- + --- J . (25) 

Finally, substituting the expansions f[2"3"j) and into Eq. (IT51) we obtain that f{u) —> 
as u — > 0, as required. 

Let us now explore the behavior of the solution as u — > oo. In this limit, Eq. (|20|) 
gives uV'{u) — >■ 0, and moreover 

2 _ nV^n) 



whence 



From fl26l) . we also find 



» v 'M~i-!S? + -' (26) 



"M~C-- + - 2 - ; + ..., (28) 



w 9m 3 

where C = const. Inserting the expansions fl26|) and fl2~Tj) into Eq. fl2TT) . we obtain 



-±^7 + 



Aq(u) u u 2 
and hence 



A (u) ~C uexp( =f- + • ■ • ) . (29) 



it 

The case C ^ is unsuitable, since Eq. f[2"9~j) would imply that A («) — > oo as u — > oo 
and, in view of formulas and (fTo"]) . the solution /(u) appears to be unbounded, 
which contradicts our assumption. Therefore, Co = and hence A (u) = 0, thus 
reducing the expansion (TT5]) to 

/(it) ~ e(Ai(«) + eA 2 {u) + • • •) exp( e - 1 y(«)). 

Arguing as above, we successively obtain A\{u) = 0, A 2 {u) = 0, etc. This indicates 
that f(u) = 0, which was our aim. 
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3. Analytical proof 

In this section, we demonstrate how the theory of g-difference equations (see Refs. 
[U [3]) can be used to show that Eq. ffl3|) has no nontrivial bounded solutions. In 
what follows, we assume that q ^ 1. 

As explained in the Introduction (see Eq. (JE])), Eq. ( Fl3l) splits into two (similar) 
one-sided equations, so it suffices to consider the boundary-value problem 

y'{x)+y(x)^ -y(qx) + ^y(q' 1 x), x > 0, 
j/(0) = 0. 

Assume that y(x) is a bounded solution of Eq. (1501) . < 5 (a; > 0), and let ?/(s) 

be the Laplace transform of y(x), 

POO 

y(s) := / e~ sx y(x)dx, 
Jo 

then y(s) is analytic in the right half-plane, 3?s > 0, and 

\y(s)\ ^>o. (31) 

On account of the boundary condition y(0) = 0, Eq. (1301) transforms into 

(1 + s)y(s) = ^((T 1 *) + (32) 

or, after the substitution y>(s) := sy(s), 

(1 + s)v(s) = ^(q- 1 s) + ^(qs). (33) 

Note that the estimate (1311) implies 

1^)1 Rs>0, (34) 

and in particular <^(s) is bounded in the vicinity of the origin. 
Let us rewrite Eq. (133]) in the form 



V (q 2 s)- 2(1 + qs) V (qs) + V (s) = 0. (35) 

Equation (1351) is a linear g-difference equation of order 2. According to the general 
theory of such equations (see Ref. [1]), the characteristic equation for Eq. (13"5T) (in the 
vicinity of s = 0) reads 

q 2p -2q p + l = 0, 
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and pi ( 2 = is its multiple root. The corresponding fundamental set of solutions to 
Eq. (1551) is given by 



<Pi(s) = Pi(s), 

mq 

where Pi{s), P2(s), and Ps(s) are generic power series convergent in some neighbor- 
hood of zero, and Ins denotes the principal branch of the logarithm. Note that the 
solution ^2(5) is unsuitable because it is unbounded near s = (see Eq. (|34|) ). On 
the other hand, the function <pi(s) is analytic in the vicinity of zero and, moreover, it 
can be analytically continued, step by step, into the whole complex plane C by means 
of Eq. ( 1551) . For instance, if q > 1 then the analytic continuation from a disk \s\ < a 
to the bigger disk \s\ < qa is furnished by the formula (p(qs) = 2(1 + s) (p(s) — <^(g _1 s) 
(see Eq. ( 155|) ). and so on. 

That is to say, ipi(s) can be extended to an entire function ip(s), which by con- 
struction satisfies Eq. ( 1551) for all s G C. But then, according to one result by Mason 
[2T] . the entire function ip(s) must be of zero order, and consequently (see Ref. (274 
§8.7.3]) it is unbounded on any ray (in particular, for s G (0, 00)), unless it is a 
constant. However, the unboundedness for real s > contradicts the estimate (1511) . 
Hence, ip(s) = const, so that y(s) = const - s' 1 , and by the uniqueness theorem for the 
Laplace transform this implies that y(x) = const, i.e, y(x) = y(0) = 0, as claimed. 

4. Probabilistic proof 

In this section, we give a probabilistic interpretation of Eq. ( 1T51) via a certain ruin 
problem, and prove that the corresponding solution is constant using elementary 
probabilistic considerations. Although our argument does not cover the whole class 
of bounded solutions, it contains some ideas that we will use in the second half of the 
paper to give a complete proof of our general result. 

Let us consider the following "double-or-half" gambling model (in continuous 
time). Suppose that a player spends his initial capital, x, at rate v per unit time, 
so that after time t he is left with capital x — vt. However, at a random time r 
(with exponential distribution), he gambles by putting the remaining capital at stake, 
whereby he can either double his money or lose half of it, both with probability 
1/2. After that, the process continues in a similar fashion, independently of the 
past history. If the capital reaches zero and then moves down to become negative, 
this is interpreted as borrowing, so the process proceeds in the same way without 
termination. In that case, gambling will either double or halve the debt, and in 
particular the capital will remain negative forever. 

More generally, if X t denotes the player's capital at time t > 0, starting with the 
initial amount Xq = x, then the random process X t moves with constant negative 
drift (— v), interrupted at random time instants o~i by random multiplication jumps 
from its location X{ = X ai -o (i.e., immediately before the jump) to either qxi or q~ 1 Xi 
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(q 7^ 1), both with probability 1/2. We assume that the jumps occur at the arrival 
times 0i,<72, ... of an auxiliary Poisson process with parameter A > 0, so that the 
waiting times until the next jump, Tj = Oi — a^i (cr : = 0), are independent identically 
distributed (i.i.d.) random variables, each with the exponential distribution 

P{r > t} = e~ xt , t > 0. 

According to this description, (X t ) is a Markov process, in that the probability law 
of its future development is completely determined by its current state, but not by 
the past history ("lack of memory") (see, e.g., Refs. [9l 110]). 

We are concerned with the ruin problem for this mo delH Namely, consider the 
probability fo(x) of becomong bankrupt starting with the initial capital x, 

f Q (x) := Pjliminf X t < \X = x\ = P X {T < oo}, x G R, 

where T := min{t > : X t < 0} is the random time to bankruptcy and P^ denotes 
the probability measure conditioned on the initial state Xo = x. 

From the definition of the process X t , it is clear that if x < then Tq = and so 
fo(x) = 1. For x > 0, we note that if the first jump does not occur prior to time x/v, 
then the process will simply drift down to 0, in which case T = x/v < oo. Otherwise 
(i.e., if a jump does happen before time x/v), the ruin problem may be reformulated 
by treating the landing point after the jump as a new starting point (thanks to the 
Markov property). More precisely, by conditioning on the first jump instant <j\ (= n) 
and using the (strong) Markov property, we obtain 

/•x/v 

fo(x) = P x {0! > x/v} + / Ae- As P X (T <oo\a 1 = s)ds 

Jo 

= e~ Xx/v + X [ XV e ~ Xs \ [fM x - vs )) + /o Or'O* - vs ))) ds 

= e~ Xx/v ^(^"(/oM + Zofrt)^ , (36) 

where in the last line we have made the substitution u = x — vs. The representation 
fl36l) implies that the function fo(x) is continuous and, moreover, (infinitely) differen- 
tiable, and by differentiation of Eq. ( |36l) with respect to x, it follows that the function 
V — fo(%) satisfies the generalized pantograph equation (cf. Eq. ( flBl ) 

v 11 

jy'(x)+y(x) = -y(qx) + ~y{q~ x). (37) 

It is easy to see that, in fact, this equation is satisfied on the whole As 
we have mentioned, fo{x) = 1 for all x < 0, and it is now our aim to show that the 

A similar ruin problem for the process with deterministic multiplication jumps of the form 
Xi i— ► qxi {q > 1), was first considered by Gaver [12] . leading to the equation y'(x)+y(x) = y{qx) (i.e., 
with advanced argument, cf. Eq. ©). The systematic theory of general processes with multiplication 
jumps was developed by Lev [TS] , 
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same is true for all x > 0, which would mean that the solution y = fo(x) to equation 
([TBI is a constant, fo(x) = l,iel 

To this end, note that the position of the process after n jumps is given by 

X ffn = {X 0n _ 1 -VT n )q in , n = l,2,..., 

where £ n 's are i.i.d. random variables taking the values ±1 with probabilities 1/2. By 
iterations (using that X — x), we obtain^ 

x an = (((z - vr i) <? ?1 - vt 2 ) q& t>r n J q in 

= (x - VT X ) _ VT2 g &+«~K„ wr n g € " (38) 

where S n := £i + £2 + • • • + £n, Sq :— 0. Note that S n can be interpreted as a (simple) 
random walk, which in our case is symmetric (i.e., P{£j = 1} = P{£j = —1} = 1/2) 
and therefore recurrent (see, e.g., Ref. [10J). In particular, the events A n := {S n -i = 
0} {n = 1,2,...) occur infinitely often, with probability 1. Furthermore, setting 
B n := {r n > 1}, we note that the events A n R B n (n — 1,2, . . .) are conditionally 
independent, given the realization of the random walk {Sk, k > 1}. Since the random 
variables r n (and therefore the events B n ) are independent of {Sk}, we have, with 
probability 1, 

00 00 

^2 P ( A n n B n I {S k }) = ^ P ( B n) U n 
n=l n=l 

— e_A #{ n : A n occurs} = 00, 

where denotes the indicator of an event. Hence, Borel-Cantelli's lemma (see, 

e.g., Ref. [10]) implies 

P(v4 n fl B n occur infinitely often | {5*^}) = 1 (a.s.), (39) 

and by taking the expectation in Eq. fl39|) (with respect to the distribution of the 
sequence {S^}), the same is true in the unconditional form, 

P(A n n B n occur infinitely often) = 1. 

As a consequence, the terms in the random series 

00 

n=l 

will infinitely often exceed the value 1, all other terms being nonnegative. Therefore, 
the series diverges to +00 (a.s.), and from (I3"5I) it follows that 

liminfX^ < (a.s.). 

n^oo 

In turn, this implies that To < 00 (a.s.), and so fo{x) = 1 for all x > 0, as claimed. 

8 Similar random sums as in Eq. psp arise in products of certain random matrices in relation to 
random walks on the group of affine transformations of the line (see Ref. [T5]V 
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5. Jump diffusions 



We now pursue a more general (and more systematic) approach. Equations of the 
form (jSJ) are linked in a natural way with certain continuous-time Markov processes 
(more specifically, diffusions with multiplication jumps). To describe this class of 
processes, let us consider a Brownian motion B^' v , starting at the origin, with diffusion 
coefficient k > and nonpositive (constant) drift — v < 0, 

dB? v = KdB t - vdt, B^ v = 0, 

or equivalently 

B? v = KB t -vt, t> 0, 

where B t = B] ,Q is a standard Brownian motion (with continuous sample paths). We 
assume that n 2 + v 2 > 0, so that B^' v does not degenerate to a (zero) constant. 

The random process B^' v determines the underlying diffusion dynamics for a pro- 
cess with jumps, (X t , t > 0), which is defined as follows. Suppose that the jump 
instants are given by the arrival times G\, Gi, ■ ■ ■ of an auxiliary Poisson process with 
parameter A > 0, so that Tj = <7j — ovi (i = 1, 2, . . . ) are i.i.d. random variables with 
exponential distribution, 

P{ r >t} = e~ xt , t > 

(we set formally <jq := 0). Furthermore, suppose that the successive jumps are deter- 
mined by the rescaling coefficients on of the form on = e&, where £j's are i.i.d. random 
variables. Then, the (right-continuous) sample paths of the process X t are defined 
inductively by 

= f x + B?' v , = (T < t < o-i, 

* ~ \ e^X CTi _ + B?' v - 0, < t < a i+1 , t = l,2,... 

That is to say, the process (Xt, t > 0) starts at point x and moves as X t = x + B*' v 
until a random time <J\, when it jumps to a random point 

X CT1 = e^X CTl _ = e 5l (^ + ) = e* l z + e §1 Ci, 

where £i := = — -B^. Thereafter, the process proceeds in a diffusive way 
as Xt = X ai + B£' v — B£ ,v until a random time 02, when it makes the next jump to 



X a2 = e^X CT2 _ = e*>{X n + B% - B%) 
= e 6+fl x + e 6+?1 Ci + e 6 C 2 , 

where £2 : = B^ — B^f, and so on. Iterating, we obtain that the n-th jump, occurring 
at time a n , lands at the point 

X Un = e«-X^_o = (X^ + B% - B^ ) (41) 
= e Sn (x + Ci + C2e" Sl + • • • + Cne" 5 - 1 ), 
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where 

£ : = B K f - B K f 

^ n cr n T n -1 (42) 

= n(B an - B a „_i) - v(a n - cr n _i) (n = 1,2, . . .) 
is a sequence of i.i.d. random variables. 



6. ^-harmonic functions 

Let us study the properties of the process X t in greater detail. 

Definition 1. The infinitesimal operator (generator) £ of a Markov random process 
(X t , t > 0) is defined by 

(£/)(,) . Ita E -' / ™'- /W , (43) 

/i — >0+ ft, 

with the domain T>(£) consisting of functions / for which the limit in Eq. (T4"3l) exists. 

In a standard way by considering possible scenarios for the process (X t ) up to an 
infinitesimal time h (see Ref. [9]), one obtains the following. 

Proposition 1. For the random diffusion with jumps (X t ) defined above, its generator 
£ acts on bounded C 2 -smooth functions as 

(£f)(x) = - vf'(x) + A(E[/(e«s)] - f(x)), (44) 

where £ is a random variable with the same distribution as any one of the i.i.d. random 
variables £i, £2, • • • • 

Definition 2. A function / is called £-harmonic if £f = 0. 

Note that the expectation in Eq. (I4"4"j) can be written as a Stieltjes integral, 



POO 

E[/(e ? x)]= / f(zx) dF(z) 
Jo 



where F(z) is the cumulative distribution function of the random variable Z = e^, 
i.e., F(z) := P{e^ < z} = P{cJ < ln^} (0 < z < 00). Hence, the equation £f = 0, 
with £ given by Eq. (T44|) . is equivalent to 

2 poo 

-yf(a) + «/» + A/(x) = A jf /(zx) dF(z), 

which is a balanced generalized pantograph equation of the form ([6]). 

Our aim is to study the class of bounded £-harmonic functions. Let us denote by 
|| • || the standard sup-norm on R: 



:= sup |/(:e) I . 



As a first step, we estimate the derivative of an £-harmonic function. 
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Proposition 2. Suppose that v > 0. If \\f\\ < oo and Cf = then \\f'\\ < oo. 
Proof. If k = then, according to Eq. (144")) . the equation £/ = takes the form 

«/'(*) = A(E[/(e«x)] -/(*)), 

which gives 

H/'ll < ^ H/ll < oo. 
For k > 0, the condition £/ = is equivalent to 

f'(x)-if(x) = -g(x) i (45) 

where 

7 :=^>0, g(x):=^(E[/(e^)]-/(x)). 
Solving equation (j4"5"]) . we obtain 

f\x) = / g( u ) e - 7( "- x) du + Ce 

J X 



■yx 

(46) 



■yx 



poo 

= / g(u + x) e~ JU du + Ce 
Jo 

where C = const. Since \\g\\ < 4Ak~ 2 ||/|| < oo, Eq. ()46l) implies that if C ^ then 

f'(x) = 0(1) + Ce rx -> oo (a; -> +oo), 

so that lima-^+oo /(x) = oo, which contradicts the assumption ||/|| < oo. Therefore, 
C = and 

11/11 < |b|| e" 7 "du = M < oo. 
The proof is complete. □ 

7. Proof of the main results 

Let (jF t , t > 0) be the natural filtration generated by the process (X t ), i.e., T t = 
a{X s , s < t} is the minimal a- algebra containing all "level" events {X s < c} (c G R, 
s <t). Intuitively, ^ is interpreted as the collection of all the information that can be 
obtained by observation of the random process (X 8 ) up to time t. As is well known 
(see, e.g., Ref. [9], Ch. 4]), if a function / is ^-harmonic then the random process 
f(X t ) is a martingale relative to (J 7 t), i-e., for any < s < t, with probability 1, 

E x [f(X t )\F s ] = f(X s ). (47) 

In words, this means that if we are trying to predict the mean value of the martingale 
f(X t ) at some future time t using its past history up to time s < t, then the best 



15 



estimate is given by f(X a ), i.e., the value of the process at the latest available time 
instant s. 

In particular, taking expectation of both sides of Eq. (1471) at s = gives 

E x [f(X t )] = fix), t > 0. (48) 

In addition, if / is bounded then, by Doob's optional stopping theorem (see, e.g., Ref. 
[221 §8]), Eq. (USD extends to 

E x [f(X T )] = f(x), (49) 

where T is a random stopping time, i.e., such that {T < t} G Tt for each t > 0. 
In words, one should be able to decide whether the random time T has occurred by 
observing the process up to a given time. The martingale property (j4"9l is crucial in 
the proof of our main theorem below, where we will apply it to the special sequence 
of stopping times. 

Theorem 3 (cf. Theorem Q}. Assume that E[f] < and P{£ ^ 0} > 0. If Cf = 

and ll/H < oo, then f(x) = const, 

Proof. For r > 0, set N r := min{n : S n < — r}, where S n = £i + • • • + £ n , and define 

oo 
n=l 

where (<Tj) are the time instants of successive jumps (see Section [3]). The assumption 
E[£] < implies (see, e.g., Ref. [10]) that liminf^oo S n = — oo (a.s.), so that N r < oo 
(a.s.) and therefore the random variable T r is well defined. 

Note that T r is a stopping time for the random process (X t ), i.e., {T r < t} G Tt 
for each t > 0. Indeed, suppose that we are given a sample path of the process X s 
up to time t, and in particular we know the time instants cij and the magnitudes of 
all thejumps prior to t. Then, using Eq. fj4T|) . we can reconstruct the corresponding 
valuesj^j = hi(X a J 'X ai - ). In turn, this allows us to determine if the threshold (— r) 
has been reached by the associated random walk S n and, therefore, whether or not 
the condition T r < t holds, as required. 

Now, applying the optional stopping theorem in the form (T4"9"|) . we obtain 

f[x) = E x [f(X Tr )] = E[f(e SNr x + terms independent of x)\ . (51) 
First, suppose that v > 0. Differentiation of Eq. flBTl) with respect to x gives 

/' (x) = E[e SN rf'(e SN r X + •■■)], (52) 
whence, using Proposition (2], we get 

ll/'ll < e^ll/'ll < oo. 

9 There is a slight problem if X ai _ = X ai = 0, but this only happens with zero probability, so 
may be ignored. 
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Letting here r — > +00, we conclude that ||/'|| = and so f = const. 
If v = then, differentiating Eq. (|52p we get 

f"{x) = E[e 2SN rf"{e s »r X + -..)], 

so that 

11/1 < e- 2 1|/l < 00. (53) 

Equation f|45l) (with 7 = 0) implies ||/"|| < \\g\\ < 00, so taking r — > +00 in Eq. fl53|) 
yields /"(^) = 0. Therefore, /(x) = c\x + Co, but since ||/|| < 00, we must have 
Ci = 0, so that f(x) = c = const. □ 

Proposition 3. Denote Aoo := jliminf^oo X t = +00}, and consider the probability 
of the event as a function of the initial point of the process X t , 

foo(x) :=P x (A oo ) = P(A oo \X = x), xeR. 

Then the function foo{x) is C-harmonic, i.e., 

{Cf oo ){x) = 0, xeR, 

where £ is the generator of the random process (X t ) given by Eq. fi44\ ). 
Proof. Conditioning on (< +00), by the Markov property we obtain 

f 00 (x) = E x [P x (A 00 \X h )} 
= E x [P Xh (Aoo)] 

whence, by the definition (j4"3~l) . it readily follows that Cfoo = 0. □ 

Theorem 4 (cf. Theorem [2]). Suppose that E[£] > 0. Then the function fao(%) is a 
nontrivial bounded C-harmonic function; in particular, foo{x) — > as x — > —00 and 
foo(x) — > 1 as x — > +00. 

Proof. The function foo{x) is ^-harmonic by Proposition [31 In order to obtain the 
limits of foo( x ) as a; ^ ±00, note that 

1 - /oo(x) = Pjliminf X t < +00) (54) 

= P x \ liminf X CTn < +00 \ 

= lim P x {X CTn < M infinitely often}. 

il/— >oo 

According to Eq. (jHJ), the condition X an < M can be rewritten as 

x + Ci + C2e" 51 + ■ • • + Cne -5 - 1 < Me~ Sn . (55) 
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Since E[£] > 0, the strong Law of Large Numbers implies that, with probability 1, 

S n ~ h.E[£] — > +00 (n — ► 00). (56) 
It follows that if the inequality (|55[) holds for infinitely many n, then 

x + liminf (Ci + C2e" 51 + • • • + Cne" 5 "" 1 ) < lim Me~ Sn = 0. (57) 

n— *oo n— >oo 

Moreover, using Eq. (|56|) and recalling that Q are i.i.d. random variables (see Eq. (|42[) ). 
it is easy to show (e.g., using Kolmogorov's "three series" theorem, see Ref. [10]) that 
the random series 

00 

n=l 

converges with probability 1. Therefore, from Eqs. (1551) . (1571) and (1581) it follows that 
for any M > 0, 

P x {X an < M infinitely often} < P{rj < -x}. 
Returning to Eq. (I5"3|) . we deduce that 

1 > foo{x) > 1 - P{v < -x} 

= P{i] > —x} — > 1 (x — > +00). 

On the other hand, writing the left-hand side of Eq. (|55p as x + r\ + 8 n , where 
<5n — > (a.s.), we have, for any e > 0, M > 0, 

P{a; + 1] < —e} < P{x + 77 + 5 n < for all n large enough} 
< P{x + r] + S n < Me~ Sn infinitely often}, 

which in view of Eq. (j54p implies 

0< foo{x) <l-P{v + x <-e} 

= P{r] > —x — e} — > (x — > —00). 

Thus, the proof is complete. □ 
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